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ABSTRACT 

Cosmic reionization holds the key to understand structure formation in the Universe, and can 
inform us about the properties of the first sources, as their star formation efficiency and escape 
fraction of ionizing photons. By combining the recent release of Planck electron scattering 
optical depth data with observations of high-redshift quasar absorption spectra, we obtain 
strong constraints on viable reionization histories. We show that inclusion of Planck data 
favors a reionization scenario with a single stellar population. The mean xhi drops from ~ 0.8 
at z = 10.6 to ~ 10“"^ at z = 5.8 and reionization is completed around 5.8 ^ z < 8.5 (2- 
tr), thus indicating a significant reduction in contributions to reionization from high redshift 
sources. We can put independent constraints on the escape fraction /esc of ionizing photons 
by incorporating the high-redshift galaxy luminosity function data into our analysis. We find 
a non-evolving /esc of rsj 10% in the redshift range z = 6 — 9. 

Key words: dark ages, reionization, first stars - intergalactic medium - cosmology: theory - 
large-scale structure of Universe. 


1 INTRODUCTION 


Cosmic reionization is one of the key events in the history of Uni¬ 
verse. Most of the available constraints on the epoch of reionization 
(EoR) come from the observations of the CMB by the Wilkinson 
Microwave Anisotropy Probe (WMAP) sate llite jKomatsu et all 
I 2 OI itiHinshaw et al.l201^:lBennett et alj20 13h. a nd from high red - 


shift (jsOs iBecker et alj200ll : IWhite et alj20^:lFan et alf2006h . 

The recent nine-year WMAP observations provide the value of in- 
tegrated Thomson sc attering optical depth Tei = 0.089 ± 0.014 
jHinshaw et alj2013l) . This in turn corresponds to an instantaneous 
reionization taking place at redshift Zreion = 10.6 ± 1.1, indicat¬ 
ing a strong need for sources of reionization at z > 10. However, 
improved measurements from three-y ear Planck mission suggest 
a lower value, Tbi = 0.066 ± 0.012, i jPlanck Collaboration et al.| 
l2015t) corresponding to Zreion = 8.8//j and therefore cuts 
dow n the demand for reionization sources beyon d redshift z = 
10 teobertson et all 1201^ : iBouwens et all l2015ah . The resulting 
reionization histories see m to explain the obser v ations of Lyman- 
Q em itters at z ~ 7 jMesinger et al.l I2OI5I : IChoudhurv et al.l 
I 2 OI 4 I) which were otherwise in tension with the WMAP con¬ 
straints. Although, these observations of cosmological data anal¬ 
ysis are based on the assumption that r eionization is a sudden 
and instantaneous process, re ce nt studies 1 BaAana^^oe^ 200lJ; 
IChoudhurv & Ferr^ I2OO6 IChoudhurvI l2009l : IPritchard et al.l 
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I2OIOI: iMitraetalJli^ I2OI2I: lOhara et al.l[ 20 l^ clearly support 
a more extended process spanning the redshift range 6 < z < 
15 . Several semi-analytical models have been proposed with a 
combination of different observations to put tighter limits on the 
reionization redshift and other quantities related to reionization 


20061; Diikstra et al. 20071; ISamui et al. 

2007 

; Iliev et alJ 

2008; 

Mitra et alj|201 ll; iKulkami & ChoudhurvI 

2011 

; iMitra et al.l 

2012; 

Cai et alj2014t). 


Another observation set that could be used to check the con¬ 
sistency of such mode ls is Fuminosity Function (UF) of high-z 


(6 < z < 10) galaxies dBouwens & Illingworthll20061; Oesch et al.l 

2 OI 2 I; [Bradley et alj 

2 OI 2 I; lOesch et alJl20l4 Bowler et alJl20l4 

McFeod et al.l 2014; 

Bouwens et alj|2015bt). This procedure has to 


deal with the yet poorly understood escape fraction of ionizing pho¬ 
tons (/esc). Despite of numerous impressive efforts in both observa¬ 
tional and theoretical studies, this quantity, as a function of galaxy 
mass and redshift, remains largely uncertain dPemandez & ShuUl 
l 201 lh . Available studies generally lead to a broad range and as 
well as contradictory trends of /esc on redshift. For example, 
iFinkelstein et al.l ( 12012h estimated ave rage /esc to be 30 % in or¬ 
der to get a fully ionized IGM at z = 6 . lKuhlen & Faucher-Giguera 
( l 2012 h found a strong redshift evolution of escape fraction in¬ 
creasing from cvj 4 % at z = 4 to unity at higher redshifts in 
order to simultaneously satisfy reionization and lower redshift 
Lyman-a forest constraints. Based on a robust statistical analy¬ 
sis on full CMB spectrum and qu asar data and usin g the obser¬ 
vations of high-z galaxy FFs, we iMitra et alj|2013l) derived that 
mean value of /esc is moderately increasing from 7% at z = 6 
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to 18% at 2 = 8. This increasing behavior of /esc on redshift 


is somewhat similar to that obtained or assumed in Inoue et al. 

2006 

1: iRazoumov & Sommer-LarsenI (l201Ch; iHaardt & Madau 

2011 

,l2012h; Ferrara & LoebI (20131); iFinlator et alJ (l2015b. More 


recently, using a hig h-resolution cosm ological zoom-in simulation 
of galaxy formation, iMa et alj ( 1201 Sh found a considerably lower 
(< 5%; much less than required by reionization models) and a 
non-evolving escape fraction. This unchanging tr end of fesc is 
as well consistent with many other earlier results (lGnedinll2008l : 
lYaiima et al.ll201lh . A decreasing tendency of escap e fraction with 
redshift has also been re ported in the literature dWood & LoebI 
l2000l(lKimm & Cenll2014l) . 

For all these reasons, here we revise our reionization models 
dMitra et al.l2013h in the light of recentiy available Planck data and 
improved measurements of high-z LFo 


2 DATA-CONSTRAINED REIONIZATION MODEL 


Let us first summarize the main features of the semi-analytical 
model for i nhomogeneous reionization u sed in this analysis, which 
is based on ichoudhurv & Ferraral d2005h and IChoudhurv & Ferraral 
d2006bl) . The model tracks the ionization and thermal evolution of 
all hydrogen and helium regions separately and self-consistently 
by adopting a lognormal probability distribu tiorfl at low densities, 
chang ing to a power-law at high densities dChoudhurv & Ferraral 
I 2 OO 5 ). The model considers the inhomogen eities in the IGM ac¬ 
cordi ng to the description given by MHR dMiralda-Escude et al.l 
|2000|) . in which once all the low-d ensity regions are ionized, reion¬ 
ization is said to be complete (see Ichoudhurvir2009ti . 

The sources of reionization are assumed to be stars (metal- 
free PopIII and normal PopII) and quasars. The contribution of 
quasars are incor porated here based o n their observed luminosity 
function at 2 < 6 dHopkins et al.ll200% and they h ave insignificant 
effects on IGM at higher red shifts (but also see rGiallongo et al.l 
I 2 OI 5 I : [Madau & HaardJl2oT5h . Furthermore, the model calculates 
the suppression of star formation in low-mass haloes (radia¬ 
tive feedback) through a Jeans mass prescription, which is com¬ 
puted self-consistently from the evolution of the thermal prop¬ 
erties of IGM throu gh the minimum circular v elocity of haloes 
that are able to cool (lQioudhur^_&Fe2^ 12005 ; also comparable 
to the simulations by Okamoto et al.l 2008 ; Sobacchi & Mesingetl 
l2013h . Our model also accounts for the chemical feedback (PopIII 
—» PopII transition) using merger-tree based genetic approach 
dSchneider et alJl2006l) .The model computes the production rate of 
ionizing photons in the IGM as 


f('ph(^) — ribfVic 


d/coll 

dt 


( 1 ) 


where, /coii is the collapsed fraction of dark matter halo, rtb is the 
total baryonic number density in the IGM and Alon, possibly a 
function of halo mass and redshift, is the number of photons en¬ 
tering the IGM per baryon in stars. However, throughout this work 


^ Throughout this Letter, we assume a flat Universe with Planck cosmo¬ 
logical parameters: klm = 0.3089, Ga = 1 — Hm, klbh^ = 0.02230, 
h = 0.6774, 0-8 = 0.8159 , n* = 0.9667 and Ip = 0.2453 
dPlanck Collaboration et alj2015h . 

^ Another commonly employed form of the probability de n sity fu nc¬ 
tion (PDF) is an updated version of iMiralda-Escude et al.l d2000l) lit 
dPawlik et alj2009ll . However we checked that, the differences in these two 
PDFs are much smaller than the errors in the data considered here and thus 
the constraints are unlikely to be affected. 


we assume Aion to be independent of halo mass. This quantity can 
be written as Nion = e./escAt-,,, where e* is the star formation 
efficiency, and AA is the specific number of photons emitted per 
baryon in stars dMitra et al.l2013h : 


2.1 MCMC-PCA constraints from Planck data 

From the above model, we obtain the redshift evolution of Nion{z) 
and other quantities by performing a detailed likelihood estima- 
tion using t he Pr i ncipal Component Analysis (PCA), following 
iMitra et ^ d201lL l2012h . We assume Nion as an arbitrary func¬ 
tion of 2 and decompose it into its principal components by con¬ 
structing the Fisher matrix from a fiducial model for Nion us¬ 
ing three different datasets: (i) m easurements of pho t oioniz ation 
rates Fpi in 2.4 ^ 2^6 from [Bolton & Haehneltl d2007h and 

1 Becker & Bolton! d2013l) 0; (ii) redshift evolution of Lyman-limit 
syste ms (LLS), dAYp/dz over a wide redshift range (0.36 < 2 < 
6 ) by ISongaila & Cowij d2010l) 0; (iii ) Thomson scattering opti¬ 
cal de pth Tel using recent Planck data dPlanck Collaboration et alJ 
l2015h . We further impose fairly model-inde pendent constra i nts on 
neutral hydrogen fraction at 2 ~ 5 — 6 from lMcGreer et ^ ]( | 2 oT^ 
as a prior to our model. We choose the fiducial model to be con¬ 
stant (no redshift dependence) as it matches the above-mentioned 
data points quite accurately. When computing the ionizing radia¬ 
tion properties, we include only a single stellar population (PopII) 
and neglect the contributions from PopIII sources (thus no chemical 
feedback in the model) since such effects will be indirectly included 
in the evolution of Nion itself. 

We then take the first 2 — 7 PCA modes (eigen-modes of 
the Fisher matrix), having the largest eigenvalues or smallest un¬ 
certainties, w hich satisfy a model-independent Akaike information 
criteria (AIC: lLiddlell2007h and finally perform the Markov Chain 
Monte Carlo (MCMC) a nalysij^ using those modes; for details, see 
lMitraetal.l ( [201lll2012l) . As the Planck collaboration has not pub¬ 
lished low-^ polarization data, we include only Planck optical depth 
data in the analysis. However, if one may wish to include the full 
CMB spectra into the analysis, one should also consider variations 
in as and Us to avoid t heir possible degeneracies with the reio niza- 
tion model parameters dPandolfi et alj201 i] : lMitra et al.ll2012h . 

The MCMC constraints on the model are shown in Fig. [T] 
The fiducial model is well inside the shaded regions for all redshift 
range. One can see that, overall our model predictions match the 
observed data points quite reasonably. We find that, all the quanti¬ 
ties are tightly constrained at 2 < 6 . This is expected as most of 
the observational information considered in this work exists only 
at those redshifts. On the other hand, a wide range of histories at 

2 > 6 is still permitted by the data. The 2-a confidence limits 
(C.L.) start to decrease at redshift 2 > 13 since the components of 


^ The Fpi measure ments of iBolton & Haehneld d2007ll and 
iBecker & BoltoiJ d2013l) somewhat depend on their choice of fiducial 
parameter set. They have provided scaling relations to convert these mea¬ 
surements to any other choice of parameters. While comparing with our 
models, we have applied appropriate scaling to those for every parameter 
set considered. The Fpi data points with errorbars in Fig. [T] reflect the 
scaled measurements. 

The reason for choosing ISongaila & Cowij i201(lh dataset s over more re¬ 
cent m easurements on mean free path of ionizing photons bv lWorseck et alj 
iioil is that the former data covers a wider redshift range, thus helping us 
to get tighter constraints on reionization parameters. 

^ All cosmological parameters are fixed at their best-fit Planck value. 
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Figure 1. MCMC results: The mean value (solid lines) and its 2-cr limits (shaded regions) for various quantities related to reionization obtained from our current 
analysis with Planck data. The fiducial model (short-dashed lines) and the model constrained using WMAP9 data (long-dashed lines) also shown for com¬ 
parison. The red points with en orbars that we have used to constrain the model are taken from th e observations of photoion ization rates Fpi jBolton & Haehnel3 
l2007t ISecker & Boltonll201.3l : top-middle panel). The redshift distributi on of LLS dA^rji/ dz jSongaila & CowidEoT^ : top-right panel) and the recent mea¬ 
surements of electron scattering optical depth Tgi from Planck mission jPlanck Collaboration et alj |2015l; b ottom -left panel). W e also show the observ ational 
limits on neutral hydrogen f r action x^\(z) (bottom -right panel) f rom various measur ements bviFan et alj (filled circle)jMcGreere^^ tlOlsI) (open 

triangl e). iTotani et alj |200^)jChomockel^l J i2013ll (open circlel. fBolton et ^ i201 lh : ISchroeder et alj j2013h fonen diamondl. [QtaeTaiTI^O^ k lOuchi et alj 
i20ld) (open sauarel. [ScheificCTeraiT]^l4) (filled square). 


the Fisher matrix are zero and there is no significant information 
from the PCA modes beyond this point. 

For comparison, we also show the model constrained by the 
WMAP9 Tel = 0. 089 zfc 0.014 value an d corresponding cosmo¬ 
logical parameters dHinshaw et al.ll2oT^ . The mean evolution of 
all the quantities for this model is almost identical to the Planck 
one at z < 6; at earlier epochs they start to differ, as expected 
from the different e.s. optical depth. However, the mean model 
for WMAP9 lies within Planck’s 2-a limits only up to 2 < 7 or 
8. We have further shown a typical admixture of PopIII contri¬ 
butions with epopiii = 0.005 (e = e»/esc) to the Wion (dotted 
black lines), which is the maximum limit that can be allowed by 
the Planck Tei data. Whereas for the best-fit WMAP9 model, we 
get epopiii = 0.014. The mean evolution of Nion{z) for Planck 
closely follows the fiducial model, suggesting that an non-evolving 
Nion is well-permitted by the current data. This is one of the key 
results from this work and deserves some more insight. 

Unlike the WMAP9 model, the smaller value of Tei from 
Planck essentially releases the need for high-redshift ionizing 
sources and allows the reionization to be completed from only 
a single stellar population (PopII). The mean evolution of pho¬ 
toionization rates Fpi shows a mild increase at 2 > 6, whereas 
the WMAP model shows a relatively higher value of Tpi at early 
epochs as it still allows the contributions coming from PopIII stars, 
which are able to produce large number of ionizing photons. A 
similar trend is also found in the evolution of LLSs. Thus the fu¬ 
ture observations on LLS around these epochs may able to fur¬ 
ther discriminate between these two models. In the Tpi plot, we 
also show the relative contributions from quasars and galaxies at 
different redshifts for the mean model by long-dashed (blue) and 
dotted (green) lines respectively. From the evolution of volume 
filling factor Qhii for HII regions, one can see that reionization 
is almost completed (Qhii 1) around 5.8 ^ z < 8.5 (2- 
a limits) for Planck model. The mean ionized fraction evolves 


rapidly in 5.8 < -z 5. 10.6, in good agreement w ith other re¬ 
cent works l lGeorge et al.l201^jRobertson et al.l2015l) . whereas the 
mean WMAP model favors a relatively gradual or extended reion¬ 
ization (spanning 5.8 < « < 13). This is reflected in the evolution 
of the neutral hydrogen fraction xhi: the mean value for Planck 
model goes from xtn 0.8 to a:Hi ~ 10~^ between 2 = 10.6 
(2 = 13 for WMAP model) and 2 = 5.8. Overall our model pre¬ 
diction for xm{z) is consistent with various observational limits 
(see the caption in the Fig.[T]for references). 


3 UV LUMINOSITY FUNCTION 

From the above data-constrained reionization mode ls, we now 
derive the high -2 LFs, following iMitra et ^ l2013l The lumi¬ 
nosity at 1500 A of a galaxy with mass M and age At (= 
tz — time elapsed between the redshift o f formation 2 ' and 
redshift of observation z) c an be written as dSarnui et '^ 120071: 
iKulkami & Choudhurvll20n]) : 

LisooiM, At) = e, MZi5oo(Af) (2) 

where, e, is the star-forming efficiency of PopII stars only as we 
will restrict ourselves to a single stellar population throughout this 
analysis. The template specific luminos ity Z 1500 is computed from 
stellar population synthesis model of ISruzual & Chariot ( l2003h 
for PopII stars having different metallicities in the range Z — 
0.0001 — 0.05. Here we have inco rpor ated the mass-meta llicity re¬ 
lation given bv IPaval et alj (l2009l) and lPaval et alj ( 1201 Ol) and take 
the appropriate model with that metallicity. The lumi nosity is then 
converted to standard absolute AB m agnitude system dOke & Gunrj 
Il983l[^lkami & Choudhurvll201lh . and finally we obtain the LF, 
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Figure 2. Left panels: Evol ution of luminosity fun ction from our model for best-fit e* (black c urve) and 2-a l imits (sha ded region) at = 6 — 10. Data (red) 
points with 2-a err ors are:lBouwens et 3 i2f)15bl) (for 2 : = 6 , 7, 8 ), combined datasets from lMcL ure et a ll j201^ and lMcLeod et al.l i2014h for z = 9 and 
improved data from lOesch et ^12014 for 2 : = 10. For completeness, we also show data from lOescI^^ iT i2013l) for 2 = 9 and lBouwens et alj fcOlSbIl for 
2 = 10 (not included in our analysis, open green points). Right panel: Redshift evolution of the escape fraction /esc with l-a errors. The 2 = 10 point shows 
the lower limit on /esc at that redshift. 


^{Mab, z), from 


${Mab,z) = 


dn 


dn dLisoo 


dMAB dLisoo dMAB 


(3) 


where dn/dLisoo is the comoving number of objects having lumi¬ 
nosity within [1/1500, ii 500 + d//i 5 oo] at a redshift 2 . This quantity 
can be calculat ed from the formation rate o f haloes using our reion¬ 
ization model jChoudhurv & Ferraral200l) . 

Now, we vary e« in eq. [^as a free parameter and match the 
observed LFs with our model predictions computed using eq. min 
Fig.mC^e/r panels), we show our results for the best-fit e, with 95% 
C.L. for redshifts 2 = 6 —10. As the observations at 2 = 10 are still 
scant, we are able to determine only an upper limit of e«. The best- 
fit e, remains almost constant (~ 4%) for all redshift range. The 
number of galaxies start to decrease at fainter magnitudes produc¬ 
ing a “knee”-like shape in the LFs which shifts towards the brighter 
ends at lower redshifts (Mab ~ —12 at 2 = 10 to Mab —14 
at 2 = 6). This is due to the radiative feedback implemented in our 
model where the star formation becomes suppressed. 

Overall, the match between data and model predictions is 
quite impressive for all redshifts considered here. The uncertainties 
(shaded regions) are larger at high -2 and at the bright end of LFs. 
For lower redshifts (2 ^ 7), although our model can match the 
fainter end of the LF accurately, it slightly over-predicts the bright 
end. This ge neral tendency has already been addressed in several 
recent works l lCai et alj20l4lDaval et al.l20l4lBowler et alj2014l 
and the references therein). In particular, Cai et al. ( 12014 ) argued 
that it can be resolved by taking the dust obscuration into ac¬ 
count, which we are neglecting here. As the dust extinction was 
insignificant at earlier times, we are getting a good match for the 
brighter end at 2 > 7. Alternatively, the surveyed volumes might 
be too small to catch the brightest, rare o bjects (for a discussion, 
see ISowler et al.ll20l^ : lOaval et alJ[2oT^ . However, it is still un¬ 
clear whether thi s discrepancy ari ses from neglecting dust or halo 
mass quenching dPeng et al.l201(]h . or it is due to a mass-dependent 
e«. Thus it would be interesting to take those effects in our model 
which we prefer to leave for future work. 


3.1 Escape fraction evolution 

As a final step, having fixed e« for different redshifts, we can de¬ 
rive limits for /esc using the reionization constraints on the evo¬ 
lution of Nion (Sec. HD with Aiy = 3200 as appropriate for the 
PopII Salpeter IMF assumed here. The uncertainties in /esc can 


Redshift best-fit e* [2-a limits] best-fit /esc [2-o'limits] 


2 = 6 

0.0368 

[0.0189,0.0743] 

0.1018 

[0.0186,0.2018] 

2 = 7 

0.0376 

[0.0194,0.0773] 

0.1215 

[0.0247,0.2424] 

2 = 8 

0.0387 

[0.0170,0.0764] 

0.1234 

[0.0309, 0.2839] 

2 = 9 

0.0390 

[0.0050, 0.0785] 

0.1202 

[0.0276,0.2838] 

2 = 10 


< 0.0455 


> 0.0996 


Table 1. Best-fit values and 2-(t limits of star-forming efficiency e* and the 
escape fraction /esc obtained from the LF matching at different redshifts 
2 = 6 — 10. At 2 = 10, we only get an upper limit of e* and a correspond¬ 
ing lower limit of /esc. 


also be calculated using the quadrature method dMitra et alJl2013[) . 
We show our resulting /esc in Table [D and the right panel of Fig. 
12 We find almost non-evolving (constant) /esc of ~ 10% (best-fit) 
in the redshift range 2 = 6 — 9. For 2 = 10, we get a lower limit 
of 10% which is a very strong constraint given the uncertainties 
present in high -2 observations. 


4 CONCLUSIONS 

Over the past few years, several numerical and analytical ap¬ 
proaches have tried to constrain reionization scenarios by using 
CMB WM AP observations and Q SOs. In particular, in our ear¬ 
lier works dMitra et alJl201lLl2012l) . we proposed a detailed semi- 
analytical modeling of hydrogen reionization using the observa¬ 
tions for photoionization rates, redshift evolution of LLS and CMB 
and succeeded to produce a good match with a variety of other rel¬ 
evant datasets. We further tested the model against the observations 
of lu minosity functions from high-redshift galaxies dMitra et al.l 
l2013h . In this Letter, w e extend those works by taking th e recently 
published Tei data from iPlanck Collaboration et al.l d2015h and vari¬ 
ous new results from the observations of galaxy LFs at 6 ^ 2 ^ 10. 

• We find that, contrary to WMAP data, a constant/non-evolving 
Wion is now allowed by the Planck data. This immediately tells us 
that, reionization with a single stellar population (PopII) or non- 
evolvin g IMF is possible, i.e. the impact of PopIII stars on reion¬ 
ization JPaardekooper et alj2013l) is likely to be negligible. 

• According to Planck data, reionization proceeds quickly from 
2 ~ 10.6 to 2 « 5.8 as the mean xyii drops from 0.8 to 10~^ 
within these epochs. We find that reionization is almost completed 
around 5.8 < a < 8.5 (2 -(t C.L.). However, the model with 
WMAP data seems to favor an extended reionization starting as 
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early as 2 « 13. Thus, the inclusion of Planck data in turn in- 
dicates that most of the reionization activit y occurs at 2 < 10 
dRobertson et alJf201.5l : lBouwens et alJ2015ah . 

• From the match between the observed high-redshift LFs and 
our model predictions, we find that the best-fit values for both e* 
and /esc remains somewhat con stant with redshifts: e« at ~ 4% 
and fesc ~ 10 % for 2 = 6 — 9 llGnedinll2008l : lYaiima et alJl201ll : 
iMa et alj20ll) . We have also obtained the tightest constraint avail¬ 
able to our knowledge on on /esc (> 10%) at 2 = 10. 

Although we have focused on high-redshift LFs, one can ap¬ 
ply the same method for the evolution in lower redshift range 
3 SS 2 5. As the dust extinction becomes significant at those 

redshifts, one has to take that into account. Moreover, the addition 
of dust and/or a mass-dependent e* may resolve the problem of 
overproducing the brighter end of LFs, as stated earlier. We defer 
these issues to future work. Also, the inclusion of full CMB datasets 
from Planck into our analysis can be helpful for ruling out some of 
the current reionization scenarios. Unfortunately, the recent Planck 
data release does not include the polarization data in their likeli¬ 
hood; instead, t hey rely on the WMAP polarizat ion likelihood at 
low multipoles dPlanck Collaboration et alj|2015ri . As most of the 
constraints at 2 > 6 related to re ionization models come from po¬ 
larization data jMitra et alj2012h . we postpone such analysis to the 
next Planck data release. 
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